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TECHNICAL NOTE 3290 


ON THE MINIMIZATION OF AIRPLANE RESPONSES 
TO RANDOM GUSTS 
By Murray Tobak 


SUMMARY 


A theoretical study is made of the motions experienced by aircraft 
in reponse to sharp-edge, harmonic, and random gusts. For the sharp- 
edge and harmonic gusts, exact responses in normal acceleration and 
pitching velocity are presented for the rectangular wing flying at Mach 
number 1.2. These are compared with approximate solutions based on 
commonly used assumptions, and the validity of each of the assumptions 
is assessed. It is determined that the use of stability derivatives in 
place of indicial functions in the equations of motion does not signifi- 
cantly impair the accuracy of solutions for transient and harmonic 
response. 

The problem of alleviating the airplane's response to random gusts 
is cast in a form amenable to treatment by the Wiener optimum filter 
theory. A derivation is given of the theoretical requirements of a 
compensating-force system that minimizes a linear combination of the air- 
plane’s mean-square normal acceleration and mean-square pitching velocity. 
Results of computations are presented which indicate the system may be 
successful in causing significant reductions of both motions. 


INTRODUCTION 


For many years, aerodynami cists have studied the motions and loads 
which aircraft experience when they encounter vertical gusts in their 
line of flight. The reason for this continued interest is readily under- 
stood when it is recalled that the stresses imposed by gusts may be among 
the most severe that an aircraft structure is required to withstand . The 
gust loading condition is therefore generally considered the crucial 
factor in determining the adequacy of a structural design. Furthermore, 
since the structural design plays a large role in dete rminin g the air- 
craft’s weight, the influence of the gust loading on the structural 
design in turn serves as a limiting factor on the max-tminn range and 
speed of the vehicle. The slightly different viewpoints expressed in 
the two reasons have inspired two fields of research, correspon ding 
roughly to the categories of analysis and synthesis. In the first, the 
necessity of including the gust loading condition in structural calcu- 
lations led to the development of analytical methods by which the motions 
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and loads imposed "by specified gusts could be predicted. Second, the 
possibility of saving weight by lightening the structure inspired efforts 
to develop aerodynamic and electromechanical devices which operate to 
offset the imposed loads caused by gusts. A brief sketch follows outlin- 
ing some of the major developments in these fields. 

The earliest work in analysis is that of Wilson (ref. 1) who fol- 
lowed the classical stability theory of Bryan (ref. 2) and Bairstow, 

Jones, and Thompson (ref. 3) . In this theory it is assumed that the 
transient behavior of the aerodynamic forces and moments following sudden 
changes in airplane motion or gust velocity can be neglected. As a 
result, the equations of motion describing the airplane’s response to 
disturbances reduce to a system of ordinary differential equations which 
are gener ally readily solved. Because of its simplicity and adaptability, 
this theory has been the one most often used in subsequent analyses. 

Later theoretical developments include the introduction of the transient 
lift functions (refs. 4 and 5) and their use in connection with the oper- 
ational calculus to describe more precisely the motions and loads caused 
by gusts of arbitrarily specified structure (refs. 6 and 7 ) • The anal- 
ysis in this case is exact within the framework of the theory used to 
derive the transient lift functions, but complete solutions generally 
involve the inversion of integral equations, which often require a for- 
bidding amount of labor. A further limitation in either of the methods 
just described is the necessity of specifying mathematically the struc- 
ture of the spacewise distribution of gust velocity through which the 
airplane flies. 'It is, therefore, not possible to handle with these 
methods the more realistic situation in which the gusts are distributed 
randomly in the airplane’s path. A significant recent development which 
overcomes this limitation concerns the introduction of concepts derived 
from the field of statistical dynamics. Here, under the assumption that 
the random process is stationary and Gaussian, it is sufficient to char- 
acterize the random distribution of gusty air by a single statistical 
quantity, the correlation function. It is then possible to predict quan- 
tities like peak load factor, or maximum accelerations on a probability 
basis (refs. 8, 9, and 10). The statistical methods still require as 
fundamental quantities, however, the responses of the airplane to step 
or sinusoidal gust inputs, so that rather than being superseded, the 
analyses described previously for deriving these quantities take on added 
importance . For successful application, these methods also require a 
depiction of the statistical nature of atmospheric turbulence, and much 
effort is currently being devoted to this task (refs. 11, 12, and 13 ). 

In the field of synthesis, a wide variety of proposals has been 
advanced for alleviating the effects of gusts. These include the use of 
spoilers (ref. 14) , gaps in the wing (ref. 15 ) , and the wing’s aeroelastic 
properties (ref. 1 6 ) to destroy or offset the increments of lift caused 
by the gust velocity. Perhaps the most successful study has been that of 
Phillips and Kraft (ref. 17) . These authors adapt the classical stability 
analysis of Wilson (ref. 1) and specify that the gust velocity disturbance 
be represented by a continuous harmonic function. The resulting mathe- 
matical simplifications enable them to study a variety of flap and 
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elevator arrangements suitable for offsetting not only the increments 
in gust lift hut pitching moment as well. Their results indicate that 
with proper selection of gearing and phase between elevator and flap 
deflections it is possible to reduce significantly the airplane's response 
in both normal acceleration and pitching velocity to the harmonic gust 
input. It r emain s for further investigations to determine whether inclu- 
sion of the transient lift functions in the equations of motion or a dif- 
ferent portrayal of the gust inputs would result in appreciably different 
conclusions . 

The present paper has two purposes, the first pertaining to analysis, 
the second to synthesis. First, in view of the widespread use of simpli- 
fying assumptions and approximations in gust response analyses, we con- 
sider the most prevalent of these in the light of more precise solutions. 
Accurate numerical results are presented for the responses in normal 
acceleration and pitching velocity to step and harmonic gust inputs for 
a rectangular wing flying at near-sonic speed. These results are com- 
pared with results containing the commonly used approximations and the 
validity of each approximation is assessed. Second, we wish to adapt to 
the gust alleviation problem a powerful synthetic development in statis- 
tical dynamics, not yet as widely used by aerodynamicists as the statis- 
tical methods of analysis previously cited. This is the optimum filter 
theory of Wiener (ref. 18) . Wiener showed that if certain statistical 
information is given about the character of message and random noise 
inputs to a linear system, one can design a filter which acts to suppress 
the noise so that a mean-square error, specified to be a measure of the 
failure of the system to follow the message, is minimized. The analogous 
terms in the gust alleviation problem are obvious: The noise corresponds 

to the random distribution of gusty air, the message to the desired path 
of the airplane, and the filter to the characteristics of a controlling 
device which operates to minimize the airplane’s response to the gusty 
air. Herein, the Wiener theory is applied to the gust alleviation prob- 
lem to derive the transfer function of a control system which minimizes 
a linear combination of the airplane’s mean-square normal acceleration 
and mean- square pitching velocity. A triangular wing flying at a high 
subsonic Mach number is used as an example in numerical calculations to 
assess the effectiveness of the control system in reducing the airplane’s 
response to a specified statistical distribution of atmospheric turbulence. 


NOTATION 


A aspect ratio 


c L 


lift coefficient. 


lift 

IcoS 


Cm 


pitching -moment coefficient. 


pitching moment 
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I 

airplane pitching moment of inertia 


I.p. 

imaginary part 


M 

Mach number 


R.P. 

real part 


S 

wing area 


V 

flight speed 


W 

gust autocorrelation function (eq. (38)) 


a 

weighting parameter in definition of mean-square error (eq. (33)) 

c 

wing root chord 


e 

base of natural logarithms 


k 

stiffness factor (eq. ( 15 )) 


i 



m 

airplane mass 


q. 

0 c 

dimensionless pitching velocity, -y 



1 2 

dynamic pressure, — p OT V 


s 

variahle of Laplace transforms 


t 

time 


w g 

gust vertical velocity 


x,y,z 

Cartesian coordinates 


*a 

distance in root-chord lengths of aerodynamic center 
of gravity 

from center 

*o 

distance in root-chord lengths of control force from 
gravity 

center of 

a 

angle of attack 


a g 

v g 

V 


3 

‘J M 2 -! 




NACA TEf 3290 


5 


7 

e2 

i 

T] 

e 

A 

M- 


f light path angle 
mean- square error 
dimensionless inertia parameter, 

dimensionless mass parameter. 


IV 2 
QboSc 3 
mV 2 


%£ c 


angle of pitch 


reduced frequency, £ x (angular frequency) 
damping factor (eq_. (15)) 


mass density of free stream 
cp chord lengths of travel in time t. 

Superscripts 



( )* complex conjugate 


Vt 

c 



mean value, defined as 


ft 

f(<p) h f 

-ft 


When a, a, and q. are used as subscripts a dimensionless derivative 
is indicated, and this derivative is evaluated as the independent variable 
approaches zero. Thus 



d ->o 
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ANALYSIS 

Airplane Motions in Response to Sharp -Edge Gust 


Consider an airplane which moves forward in level steady flight and 
suddenly encounters a uniform field of vertical velocity. The "boundary 
of the field is oriented parallel to the airplane* s lateral axis and 
extends on either side "beyond the widest span of the airplane. This is, 
of course, the mathematical idealization of the sharp-edge gust. As the 
airplane penetrates the gust front, loading develops on that portion of 
its surface influenced by the gust. The loading may be resolved into 
time-varying force and moment excitations which disturb the airplane’s 
equilibrium and cause it to undergo a plunging and pitching motion. Our 
purpose will be to calculate this motion as a function of chord lengths 
of travel cp, letting cp = 0 correspond to the instant when the foremost 
point of the airplane just penetrates the gust front. 

Coordinate system .- In sketch (a) is shown the position of the air- 
plane at an instant when only a portion of the airplane has penetrated 



the gust front. As shown in the sketch, an x,y,z coordinate system is 
placed in the airplane, the origin being fixed at the center of gravity. 

The instantaneous direction of the path of the center of gravity relative 
to still air is indicated by the direction of the flight velocity vec- 
tor V; we assume that the magnitude of V remains constant throughout 
the motion. The positive branch of the x axis is rearward and alined 
with the chord line of the wing. The y axis is perpendicular to the 
vertical plane of symmetry and coincides with the lateral axis running 
through the center of gravity; the z axis lies in the vertical plane of 
symmetry, perpendicular to the xy plane. We define the angle of attack 
a to be the angle between the xy plane and the plane containing the 
velocity vector; the angle of pitch 9 is the angle between the xy plane 
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and the horizontal plane. The flight path angle 7 is the angle "between 
the horizontal and the plane contain in g the velocity vector; it is equal 
to 0-a. Finally, the gust velocity is indicated "by w g , and is measured 
relative to still air. All quantities are indicated in sketch (a) in 
their positive directions. Forces are measured positive upward, whereas 
pitching moments are positive when tending to increase the angle of pitch. 

Equations of motion .- In this analysis, the indicia! function concept 
will be used to write the equations which define the airplane motion. 

Since this will involve the use of both the superposition principle and 
the results of linearized unsteady- flow theory, the applicability of the 
analysis will be limited by the conditions under which their use can be 
justified. Hence, for example, disturbances and angular displacements 
must be small. It follows that the equations of motion may be written 
as 


mV 2 


7* (cp) = 2 forces 


Of 


q’ (cp) =2 pitching moments 


( 1 ) 


where the independent variable is chord lengths of travel cp. The forces 
and moments to be summed are of three types, those caused (l) by varia- 
tions of the angle of attack, a(cp); ( 2 ) by variations of pitching veloc- 
ity, q.(cp) j and (3) "by the gust velocity, w g . If the indicial lift and 
moment responses corresponding to step changes in each of the variables 
a(cp) , q(cp) , and w g are given, the forces and moments due to variations 
in these quantities may be built up by use of the superposition integral 
(cf., for example, ref. 19) . For the case w g = constant, equations (l) 
may be written as 


T)7’ (cp) 

Sq.* (<p) 


_d 

dcp 


cp cp 

J Cr tt (S)a(q>-6)di + ^ / ^q^)a(<p-i) d i + ^ g C Lg (cp) 


> ( 2 ) 


- ~ J P-6)di + ^ / C mq(£)<l(cp-0<li + l^gGmgCcp) 


where 


CLaM, W?) 

CL a (cp)^ Qmq(<P) 


indicial lift and pitching-moment responses to • uni t step 
change in angle of attack, a 

indicial lift and pitching-moment responses to unit step 
change in dimensionless pitching velocity, q = kh 
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CL g (<P)> Cm g (<P) 
lag 


indicial lift and pitching-moment responses to unit step 
change in dimensionless gust velocity 


step change in dim ensionless gust velocity. 


Is 

V 


Equations (2) are a pair of integral equations involving the unknown 
variations a(cp) , q.(cp) , y x (cp) . Any one of these may "be eliminated by use 
of the equality 0 = a+7. Vie shall "be interested mainly in the normal 
accelerations and pitching velocities experienced by the airplane, and 
hence choose to eliminate a. Further, it will be found advantageous to 
separate out the steady-state values of the indicia! functions. Thus, let 


Cpjcp) = Ci a (co) - Fi(cp) 

CLq(q>) = Cp^Coo) - F 2 (q>) 

C^Ccp) = Cjn^oo) - F 3 (cp) 

Qm^T) = QniqCoo) - F 4 (cp) 


*\ 


} ( 3 ) 


The quantities Cp^Xoo) f Cp (00) , 0^(00) } and 0^(00) are the steady-state 

values of their respective indicial variations, and hence are also equal 
to the stability derivatives bearing the same subscripts. Substituting 
the expressions (3) in (2), carrying the derivatives through the integrals, 
and using the facts that a(0) = q(o) = 0, a* (cp) = q(cp) - 7' (9), we get 
for the equations of motion 


T 17 * (cp) 



[q.(i)- 7 *( 6 )] 4 l + 


CL a (°°)q.(<p) - 




cp cp 

J Fx(6)[qL(9-l)-7 , (<p-l)]d| - J F 2 (|)q*(<P-l)d| + la g C Lg 


(cp) 


<P 

£q.’ (cp) = 0^(00) J [q(i)- 7 * ( 6 ) ]d| + 0^(00) q(q>) - 


> W ' 


cp 9 

J Fs(i) [qX^-SW’ (9-1) ]d| - J F 4 (£)q* (<p-&)dS + ^^(t) 


J 
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Transformed equations of motion .- Equations (4) can "be solved for 
the unkn own quantities 7* (cp) and q.(q>) hy the use of Laplace transforms. 
With the notation 


00 

L[F(cp)]= f e“ sc %-(cp)dcp 


(5) 


let 


LCr* (<p) 1 = sr(s) 

L[q.(<p)] = Q(s) 

L[PL g (<P)] = c Lg( s ) 

L [pmg(*?) ] = c mg( s ) 

L[Fi(<p)] = fi(s), i = 1,2,3,4 




> 


J 


( 6 ) 


Carrying out the transformations (6) in equation (4) and solving for sr(s) 
and Q(s), we get 


r(s). 8 

c Lg( s ) D ( s )“ c mg( s ) B ( s )] 

ln, g 

A(s)D(s)-B(s)C(s) 

Q(s) _ 8 

G m K ( s ) A ( s )“ c Lg( B ) c ( s ) ] 


A(s)D(s) -B(s)C(s) 


where 


A(s) = s[r]-fi(s)] + (^(00) 

B(s) = s2f 2 (s) - s[Cl^(co) -fi(s) ] - Ci^(co) 
c(s) = -sf 3 (s) + Cn^Cw) 

D(s) = s 2 [^+f 4 (s)] - s[C HqL (co)-f 3 (s)] - Cm^oo) 
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Stability derivative analysis .- Since we shall be interested in com- 
paring the transient responses that result from solution of the exact 
equations ( 7 ) with those resulting from approximate formulations, it is 
instructive at this point to obtain the transformed equations which derive 
from the approximate but more familiar stability derivative analysis. In 
the latter case the equations of motion (in terms of chord lengths of 
travel) may be written 

47‘(<P) = \(“) a ('P) + PLjOt* (cp) + C^co) q(cp) + lXgC Lg (cp) 

Cq.’ (q>) = Qm a (“) a -(q>) + (cp) + 0 ^( 00 ) q(cp) + ln, g Cm g (cp) 

where Cp^ and are the (constant) stability derivatives due 

form vertical acceleration. Usually, in stability analyses, the 
OLq^oo) and Cp^ are discarded; they will be retained here, however, in 

order to demonstrate fully the degree of correspondence that exists 
between equations (4) and (8). Further, in some gust analyses, the indi- 
cial gust functions Cp, (cp) and Cm (cp) are replaced by constants, equal 

o o 

to their respective steady-state values, Cp^Joo) and Cn^C 00 ) • The effect 
of this approximation will be considered subsequently. 


* (8) 

to uni- 
texms 


Taking Laplace transforms in equation (8), eliminating a as before, 
and solving for sFx(s) and Qi(s), we get 


sri ( s ) s [ c L g ( s ) D i( s )- c m g ( s ) B i( 8 )] 

In,g A 1 (s)Di(s)-B 1 (s)C i (s) 

q 1 ( s ) _ s [ c m g ( 8 )Ai(3)-c Lg (8)C 1 (s)] 
lo-g Ax(s)Dx(s) -Bx(s)Cx(s) 


(9) 


where 

Ax(s) = s(q+Cp^) + Cp^oo) 

Bx(s) = -s[Cp,^(co)+Cp^] - Cp^oo) 

Cx(s) = Qm ct ( 00 ) + s Qmd, 

b i( b ) — s 2 £ - s [Cmq(°°)+Cm^] - Cm a (co) 

and the subscript 1 is intended to distinguish the quantities from 
their exact counterparts in equation ( 7 ) . 
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In order to make clear the approximations that are implicit in the 
use of equations ( 9 ) let us now derive equations ( 9 ) from the exact equa- 
tions ( 7 ). To do this, expand the exponential e -sc P in a power series 
in s in each of the functions fi(s) in equations ( 7 ). Thus, for 
example, for the function f 3 (s) we have 


CO 00 00 

fs(s) = yV 3 (cp)dcp - s J cfF 3 (q>)dq) + s 2 J 2~ F 3 (cp)dcp 


( 10 ) 


Wow it is known (cf. ref. 20) that the first term in this expansion is 
equal to -Cj^, the pitching-moment coefficient proportional to the 
motion a* (cp) = constant. By a slight extension of the development of 
reference 20, one can go on to show that the remaining coefficients in 
equation (10) are likewise stability derivatives corresponding in order 
to the motions a M (cp) = constant, a* *’(cp) = constant, etc. The same is 
true for the other indicial function transforms, and in fact the following 
identities can be shown to hold: 




Sc T 


^ C L 


(n) 




5Qm 

aJ n) 


dq 


(n) 


= (-l) n f 7^7 Fi(q>)dcp } n>l 
J 0 (n- 1 ) ! 

= (-l) n [ -r~— F a (cp)dcp 
Jo (n- 1 ) I 

OO 

= (-l) n f F 3 (cp) dtp 

J Q (n- 1 ) ! 

= (-l) n f 7^7 F 4 (cp)dcp 
Jo ( n_ T) ’• 


> (H) 


J 


When the forms (10) and (ll) are 
of the result with equations ( 9 ) 
the result of retain i ng only the 
fi(s) and f 3 (s) . 


substituted in equations (7)* comparison 
reveals that the latter equations are 
constant terms of the expansions of 


One is. tempted to propose that the transient responses resulting 
from solution of equations ( 9 ) would agree more closely with those of 
equations ( 7 ) if some of the higher-order terms were to be included in 
equations ( 9 ) . Such a proposal is incorrect, however, for the following 
reason: Retention of terms of order s and higher raises the degree of 

the denominators of equations ( 9 ) and hence the number of roots. The 
number of requirements for stability is therefore also raised, and it can 
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le shown that these added requirements are spurious. In fact, the only 
stability requirements that are consistent with those of the exact equa- 
tions (7) are those which exist for the foims (9) • Therefore, we shall 
consider that the correct form of the stability derivative equations iB 
that which stands in equations (9) , and we shall make our subsequent 
comparisons on this basis. 1 

Solutions for 7’ (9) and q(cp).- The inversions of equations (7) give 
the complete histories 7’ (cp) and q(cp) subsequent to the airplane’s first 
penetration of the sharp-edge gust. Unfortunately, for most cases of 
interest, the equations for the indicial functions are so complex as to 
preclude attempts at direct inversions of equations (7) • Several alter- 
native techniques have been proposed. One such technique, which appears 
attractive for application to supersonic speeds, should be specifically 
warned against. Since the supersonic indicial functions reach constant 
values in a finite time, all approximations which take this behavior into 
account result in transcendental Laplace transforms of the type f(s)e" Sc; , 
where a is a constant. One is then tempted to expand the exponential 
function in a series of polynomials in s. It has been demonstrated 
(cf. refs. 21 and 22) that expansions of this sort, in which higher-order 
terms are retained, lead to fallacious results. An expansion technique 
that is valid and applicable to both subsonic and supersonic speeds is 
to put equations (7) into the form g(s)/l+k(s), expand in the series, 
g(s) [l-k(s)+k 2 (s) - . . .] and invert term by term. This gives a series 
of convolution integrals which can be shorn to be equivalent to the 
Liouville expansion used by Lomax in reference 7- The method is easy to 
apply but becomes excessively tedious if more than three or four terms 
of the series are required to insure adequate convergence. Rather than 
adopt this method, we shall make use of one recently presented by Huss 
and Donegan in reference 23- Here it is necessary to have the response 
to harmonic inputs of the system with the desired transient response. 
Having this information, one can, by comparatively simple and rapid numer- 
ical operations, extract the transient response to any desired accuracy. 
Therefore, rather than attempt to calculate the transient responses 
directly from equations (7)., we turn instead to the calculation of the 
harmonic responses. 


Airplane Motions in Response to Harmonic Gusts 


It is well known (cf., for example, ref. 24) that the harmonic 
response and the transient responses to step or impulsive inputs are 
closely related. Having given the Laplace transforms of the step responses 
in equations (7) we may write down the harmonic responses directly from 
them by means of the relations 


^Tf desired, one can retain the terms dCp/dq’ , BCm/Sq’ in equa- 
tions (9) since their retention does not raise the degree of the denomi- 
nators. Their effect, however, is negligible. 
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> ( 12 ) 



where A is the ( dim ensionless) frequency of the harmonic gust input. 


Exact and approximate solutions .- Equations (12) are valid under 
the restriction that the Laplace transforms in equations (7) sr(s)/la g 
and Q(s)/la g represent stable transfer functions. Assuming this is 
true, we may write the exact harmonic responses as 


z! 

a g 


(iA) 


-A 2 [c Ls (iA)D(iA)-c mg (iA)B(iA) 
A( i A) D ( i A) -B ( i A) C ( i A) 



Cmg(iA)A(iA)-c Lg (iA)C(iA) 
A( i A) D ( i A) -B ( i A) C ( i A) 


> 


J 


(13) 


where 

A(iA) = iA[n-fi(iA)] + 0^(00) 

B(iA) = -A 2 f 2 (iA) - iA[CL q (”)"fi(iA)] - 
C(iA) = -iAfs ( iA) + Qm a ( 00 ) 

D(iA) = -A 2 [$+f 4 (iA)] - i A [ C m ^( 00) -f 3 ( i A) ] - 


Likewise, of course, the harmonic responses as derived from the stability 
derivative analysis carry over from equations (9) : 


^ («0 


a., 


'6 




c Lg (iA)D!(iA) -c mg (iA)B 1 (iA) 
A 1 (iA)D 1 (iA) -B 1 (iA)C 1 (iA) 






-A 2 


c mg (iA)A 1 (iA)-c L (iA)Ci(iA) 
Ai(iA)D 1 (iA)-B 1 (iA)Ci(iA) 


> (1^) 


J 
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where 


Ai(lA) = IAt ) 1 + Ci^Coo) ; V = r\ + Oj^ 
Bi(iA) = -iA[CL^(oo)+Cl^] - Cl^co) 

Gi(iA) = iACj^ + 0^(00) 

Bx(lA) = -A 2 £ - iA[C ] ]^(oo)+C m( ^] - C ma (i») 


and, as before, the subscript 1 distinguishes the quantities from their 
exact counterparts in equations (13) . 

Numerical solutions .- With the use of equations ( 13 ) , it is possible 
to compute exactly the harmonic responses in normal acceleration and 
pitching velocity of those wings whose judicial lift and moment responses 
to step changes in gust velocity, angle of attack, and pitching velocity 
have been calculated. For supersonic speeds, the necessary solutions are 
available for the two-dimensional wing (refs. 25 and 26), the rectangular 
wing (refs. 27 and 28 ), and the supersonic-edged triangular wing (refs. 29 
and 36) . For use in equations ( 13 ) these solutions must be transformed 
to functions of A by means of equation ( 5 ) (with s=iA) . In order to 
facilitate the use of equations (13) , the transformations have been 
carried out for the three classes of wings . The results are compiled in 
Appendix A. 


Comparison of Exact and Approximate Step and 
Harmonic Gust Responses 


Having solutions for the exact harmonic responses and a convenient 
method of extracting from them the transient responses (ref. 23) , we are 
now in a position to study the effects of some of the approximations and 
assumptions commonly made in gust analyses. In order to examine these 
effects under conditions where they might be expected to be important, 
numerical calcul at ions for the step and harmonic gust responses have been 
carried out for a rectangular wing of aspect ratio 3 flying at Mach num- 
ber 1 . 2 . The reason for this choice of Mach number was that it is prin- 
cipally at speeds near the sonic speed that the indicial functions show 
large variations. Hence, approximations to the indicial functions would 
be expected to be least valid there. The rectangular wing was chosen for 
study because of the necessity of having available a complete set of theo- 
retical indicial functions applicable to Mach numbers near unity. It is 
recognized that for tailless airplanes the delta plan form would have 
been a more appropriate choice than the rectangular wing that was selected. 
The numerical results are therefore somewhat unrealistic; nevertheless > 
the conclusions to be dram here regarding the validity of approximate 
solutions may reasonably be expected to apply as well to other more 
practical plan forms. 
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Harmonic responses .- Figures 1 and 2 show the wing’s normal acceler- 
ation and pitching velocity responses to harmonic gust inputs as calcu- 
lated from the exact transfer functions, equations (13), and from the 
approximate stability derivative formulation, equations (l4) . The con- 
stants which define the airplane’s inertial properties are tj’ = 100, 

£ = 100; the center of gravity is located at the wing leading edge. It 
is apparent that except in the vicinity of the peaks of the curves, the 
approximate stability derivative results provide excellent representations 
of the true variations. It will be recalled that there were essentially 
two assumptions involved in the use of the stability derivative expres- 
sions; first, that the contributions of the functions f 2 (iA), f^iA) 
could be neglected, and second, that the contributions fi(iA), f 3 (iA) 
could be replaced by the constants -Cq^ and It has been determined 

that such differences as do appear between the exact and approximate 
results are caused almost exclusively by the second of these assumptions; 
hence, neglect of the contributions of the indieial functions due to 
pitching velocity is found to be .justified. As for the discrepancies 
caused by the second assumption, it should be noted that the frequency 
band over which they occur is very narrow (approximately 0.1<A<0.2). 

The differences in areas beneath the exact and approximate curves in this 
band are quite small, and hence one may anticipate that the discrepancies 
evident on figures 1 and 2 will cause only insignificant differences in 
the corresponding transient responses. 

Step responses .- Figures 3 and 4 show the transient responses in 
normal acceleration and pitching velocity to a step change in gust veloc- 
ity as calculated from the harmonic responses and the use of reference 23. 
It is clear that in the important features of the curves, namely, the 
maximum magnitudes and the points at which they occur, the differences 
between the exact and approximate results are not significant: The values 

of cp at which the maxi mum values of normal acceleration and pitching 
velocity occur are given by the approximate results almost without error; 
the error in maximum normal acceleration is less than 10 percent, whereas 
the error in maximum pitching velocity is hardly measurable. 

Approximations to gust functions .- It may be instructive to point 
out that the conclusion to be drawn from the preceding section, that the 
use of stability derivatives in place of indieial functions is valid, is 
more general than is first apparent. The denominators of equations (7) 
form the "characteristic" part of the transfer functions - that is, the 
part that does not change with the type of excitation function. Having 
show that stability derivatives may be used in the denominators without 
serious damage to the responses, we may expect that similar good results 
will be found for the responses to many other types of excitation. The 
effects of approximations to the excitations themselves, however, (in 
this case the gust functions) still remain to be studied. 

It has been customary in some gust analyses to replace the gust 
functions Cq^cp) and C m ^(q)) with step functions, equal in magnitude to 

their respective steady-state values Cq^co) and 0^(00) . Thus, in the 
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transfer functions, equations (9)> the 

and Cjj^s) are replaced "by the terms 

tage is, of course, that now equations (9) are easy to invert. The 
results for the transient responses are 


transformed gust functions c T (s) 

J-icr 


■6 


Cj^(oo) /s and 0^(00 ) /s. The advan 


where 


~ (?) = ^(cojE’ (<p) + AiR(cp) 

S 

~ (?) = 7^R(?) 


■\ 




J 


R (?) 

"n* 


^ e~^sin kcp 


V- = 


1 W* (^Qmq+Qny^ +c L q ^ c it\x- c m <1 c I^ 


ri* 




k = 


- tiCh^+Gl^Ch^- Cm^Ci^ 


-H 2 ' 


Ai = Cn^ ^CLg+Ci^ - ^Qmq+Pm^ 

^2 = li’Cma " 


(15) 


These results sire compared with the exact transient responses also in 
figures 3 end It is clear that in this case the errors, at least in 
the normal acceleration response, may he significant; the maximum value 
of normal acceleration, for example. Is overestimated hy 25 percent. The 
errors in pitching velocity (fig. 4) are not as large; however the nature 
of the response for sma l 1 values of cp is not correctly preserved, so 
that large errors will appear in the pitching acceleration response. As 
a saving feature, it is noted that the errors caused hy the substitution 
of constants for the gust functions will diminish in severity with increas- 
ing supersonic Mach numbers, since the gust function variations do in fact 
rapidly approach steps for larger Mach numbers. 
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Neglect of pitching .- Since the loads caused "by normal accelerations 
are in general more severe than those caused ty pitching, many investiga- 
tors have neglected the airplane’s pitching degree of freedom altogether 
and assumed that the airplane is free to move only in the vertical direc- 
tion. In order to investigate the validity of the foregoing assumption, 
numerical calculations have been carried out for the same operating con- 
ditions as used above, considering that the airplane is now restrained 
from pitching. 1 


The exact harmonic response may be derived from equations (13) simply 
by letting the inertia parameter £ approach infinity. The result is 



-A 2 c Lg (iA) 
iA[T]-fi(iA) ]+CiJ«>) 


(16) 


Likewise, the harmonic response as obtained from the stability derivative 
formulation (eqs. (l4)) is 


% to) - 


iA 


-A 2 ei, s (lA) 


(17) 


The corresponding transient responses may again be extracted from these 
results by use of reference 23. Numerical results for the harmonic and 
transient responses are shown in figures 5 and 6. 


Note first in both figures 5 and 6 that, as in the previous cases, 
the differences between the exact results and those obtained from the 
stability derivative formulation are not significant. Next, compare the 
transient response for the single degree of freedom motion (fig. 6) with 
the exact response when pitching is included (fig. 3) • It is noted that 
neglect of the pitching degree of freedom causes the maxi mum value of 
normal acceleration to be overestimated, but by only about 12 percent. 
Hence, use of either of the single degree of freedom equation (l6) or 
(17) to predict maximum normal acceleration is probably justified. It 
may be necessary to point out, however, that even though the effect of 
the pitching degree of freedom on the maximum value of normal acceleration 
is small, the use of equation (l6) or (17) rather than (13) or (l4) will 
generally prevent one from obtaining the true nature of the remainder of 
the transient response curve. The reason for this can be most clearly 
seen by comparing the denominator or characteristic part of equation (17) 
with that of the equation in which the pitching degree of freedom is 
included, equation (l4) . The denominator of equation (17) is of first 
degree in A, whereas that of equation (l4) '(•rill be found to be of second 
degree. The transient response corresponding to equation (17) will there- 
fore always be a subsidence whereas the transient response corresponding 
to equation (l4) will most generally be oscillatory. 
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Concluding remarks .- The major conclusion to "be drawn from the 
results of this study is that the use of the stability derivative formu- 
lation of the equations of motion does not significantly impair the accu- 
racy of the solutions for transient and harmonic responses. It should 
he pointed out that the test case chosen for study represents conditions 
under which the differences between the exact and approximate formulations 
might be expected to be most apparent. In view of the excellent agreement 
between the two results , it should be anticipated that the stability deriv- 
ative formulation will yield accurate results over a wide range of operat- 
ing conditions and for many other types of plan forms. 


SYNTHESIS 


It is clear that since most gust disturbances are random in nature , 
this fact should be considered in a realistic analysis of the airplane’s 
response to gusts. In recognition, a number of authors have studied the 
problem and have shown how the step and sinusoidal gust responses can be 
adapted to serve in the analysis of random motions as well (refs. 8, 9 , 
and 10) . The methods are based on concepts derived initially in the 
field of communications theory; their use is valid under the assumption 
that the disturbance mechanism can be described statistically aB a sta- 
tionary random process (ref. 2k ) . As applied to the gust problem, the 
assumption appears to be a reasonable one in the light of evidence avail- 
able from experimental studies of atmospheric turbulence (refs. 11, 12, 
and 13 ) . 

There is still another aspect of communications theory, however, 
valid within the same assumption and building on the same methods, which 
has not yet received wide attention in aerodynamic applications. This 
is the optimum filter theory of Wiener (ref. 18) . The theory differs in 
emphasis from those mentioned earlier, in that the aim is to synthesize 
a system which m ini mizes responses to unwanted random disturbances; in 
contrast, the aim above is to analyze the responses of a given system. 

The connection between the aim of the Wiener theory and of gust allevia- 
tion is evident. In view of the gains possibly to be realized by devising 
a gust alleviation system which takes into explicit account the actual 
random nature of gust disturbances, consideration is given below to the 
adaptation of the Wiener theory to this end. 


General Considerations 


In the following general discussion of the gust alleviation problem 
several assumptions will be made regarding the availability of the quan- 
tities which are needed in order to define and solve the problem. It is 
believed that these assumptions are reasonable; however, in view of the 
exploratory nature of the study, questions concerning practical realiza- 
bility of the quantities will not be considered. 
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Gust al 1 eviation problem .- The problem to he considered can he 
illustrated by means of the block diagram in sketch (b) . As indicated 



in the sketch, it is necessary that the airplane be equipped with a gust 
sensing device which sends to the control system a signal proportional 
to the gust velocity. It will be assumed that in response to commands 
from the control system there can be supplied a force and a moment which 
operate independently. It will also be assumed that the statistical 
nature of the gust disturbance is known, and that the responses of the 
airplane to step changes in the gust velocity and the control force and 
moment are known. Given this information, the problem is how to design 
the block labeled "control system" so that the force and moment it con- 
trols will counteract as nearly as possible the forces imposed by the 
gust so as to minimize the deviations of the airplane path from a desired 
path. 


Here, it will be desired that the airplane fly a straight path, 
without pitching. It will be considered that the airplane motions most 
pertinent to the problem are the variations in normal acceleration and 
pitching velocity. With these specifications, a more explicit represen- 
tation of the problem can "be constructed as in sketch (c). 

The block representing the control system in sketch (b) is now indi- 
cated by the two transfer functions Kp and IQ 4 which send commands, respec- 
tively, to the control force and control moment. These commands are in 
response to signals from the gust sensing device. Similarly, the block 
representing the airplane in sketch (b) is indicated by the six transfer 
functions, relating inputs in control force, control moment, and gust 
velocity to the outputs in normal acceleration and pitching velocity. 

The net airplane motions 7 ’ and q. are of course the sums of the sepa- 
rate responses to the force, moment, and gust inputs; this is indicated 
in sketch (c) by the junctions marked with positive signs. 
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Sketch (c) 

Equations of motion .- If we say that the "blocks in sketch (c) are 
representative of the responses to step changes in the various inputs, 
then for arbitrary variations of the inputs, the system’s equations of 
motion can "be written "by a straightforward application of the superposi- 
tion integral. Thus, as functions of chord lengths of travel cp, the 
airplane motions 7 ’ (9) , q(cp) are 

cp cp 

7 ’ (cp) = JL P — (cp-|)F(g) d i + f — 

d( PJ 0 IP d 9j 0 1M 

9 

[■—- (<P-6)«g(i)d| 

d 9 Jo tog 

9 9 

q(cp) =if 4 (9-l)F(|)d| + JL P JL 

dcpJo IP dcpJo 124 

9 

JL P (9-l) a g(O d i 
d 9 J 0 tog 


(cp-|)M(|)d| + 

> (18) 

(cp-i)M(i)d| + 
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Similarly, the control force and moment variations can he written 


F(cp) 

M(cp) 


9 


%(<p-6)ag(6)a6 


9 


(p-Dag(i)ai 




) (19) 




The quantities -2— (cp), -0- (<p) in equations (18) are the responses in 
■^g ^g 

normal acceleration and pitching velocity to a step change in dimension- 
less gust velocity cc,g, and hence are identical to the responses consid- 
ered in detail earlier in this report. Likewise, the quantities 2L (cp), 

TT* (<P)> A (9). “7 (9) are the responses to step changes in force and 
JLb 1M 

moment inputs; their derivation is entirely analogous to that of the gust 
responses. Finally, the quantities Kp(cp), K^(cp) in equations (19) are 
the responses in control forte and control moment, respectively, to a 
step change in the gust velocity. These quantities are the ones which 
are to he constructed such that the responses to random gusts are 
minimi zed. 


Taking Laplace transforms in equations ( 18 ) and (19), we get 




sr(s) 


s 2 r(s) 


IF 


■_ f(s) + par] ” (s) + 


2 r(s) 


la,, 


g 


V s ) 


Q(s) 

f(s) 


sQ(s)"| w % f sQ(s) ] 

if— J f(B) + ps-J m(s) + 


Sk F (s) Arr(s) 


'sQ(s) 


la, 


g 


V s ) 


Ag(: 


> ( 20 ) 


m(s) = 


3 %(s) 


Ag(B) 
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■where 


f(s) = L[F(<p)] 
m(s) = L[M(cp)] 

Ag(s) = L [cLg(cp) ] 

ki(s) = L [Kj_(cp) ] , i = F,M 


The bracketed quantities in equations (20) are the "system functions" of 
the various responses; they are equal to the transforms of the responses 
to impulsive inputs. Thus, for example, if Q(s)/lF is the Laplace 
transform of the response in q(cp) to a step change in F, then sQ(s)/lF 
is the transform of the response in q(cp) to an impulsive change in F. 

In order to simplify the notation, let the system functions he 


oi(s) , ihh E l(s) 
Ms) - — Ms) 

03(8) BsCs) 

lx g 

X(s) = skg,(s) 
Y(s) = sk M (s) 


sQ(s) 

IF 

sQ(s) 

1M- 

bQ(s) 

ICCg 




> ( 21 ) 


J 


Substituting equations (21) in (20) and solving for sr(s) and Q(s) , we 
get 


sr(a) = [Gi(s)X(s)+G- 2 (s)Y(s)+G3(s) ]Ag(s) 
Q(s) = [Hi(s)X(s)+H 2 (s)Y(s)+Hs(s) ]Ag(s) 


) ( 22 ) 


Perfect Alleviation 


Before considering the application of Wiener theory to the gust 
alleviation problem it is instructive first to study the possibilities 
of perfect alleviation. By perfect alleviation we mean that the control 
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system operates so as to reduce the airplane responses in normal acceler- 
ation and pitc hin g velocity identically to zero for a ll, includ i ng random, 
gust inputs. 

Reduction' of motions to zero.- We note in equations (22) that, 
irrespective of the type of gust input, in principle it is possible to 
reduce identically to zero the airplane’s response in both no r mal accel- 
eration and pitching velocity. Thus, for any Ag(s), sr(s) and Q(s) are 
zero if 


_ Ha(s)G a (s)-H 2 (s)G 3 (s) 
Gx(s)H 2 (s) -G 2 (s)1x(s) 

y( ^ Hx(s)G3(s)-Ba(s)Gx(s) 
" Gx(s)H 2 (s)-G 2 (s)Hx(s) 


"v 


) ( 23 ) 


On reflection, this result is obvious; in effect it simply specifies that 
the control force and moment supplied in response to a step gust must 
exactly balance the force and moment imposed by the gust. Clearly, if 
the control system is capable of canceling the step gust input it will 
cancel all other gust inputs as well, since these may be viewed merely 
as successions of steps. 

Alleviation with one control .- In view of the obvious complexity of 
a control system that might be built according to equations (23) we con- 
sider next the possibility of perfect alleviation when only one control 
is provided. Thus , , we eliminate the independent control pitching moment 
Y(s) in equations (22) and consider that only the control force (and its 
accompanying pitching moment) is available for counteracting the forces 
and moments imposed by the gust. 

With Y(s) = 0, equations (22) become 


sr(s) = [Gx(s)X(s)+G 3 (s)]Ag(s) 
Q(s) = [Hx(s)X(s)+BV 3 (s)]A g (s) 


(2k) 


Perfect alleviation is still a possibility in the event that 
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In general, this will not he the case, as can he seen if it is recalled 
that in order to provide for perfect alleviation the force and moment 
imposed by a step gust must he canceled exactly. Now the force due to 
a step gust will in general huild up monotonically from zero to a steady- 
state value whereas the pitching moment will change sign some time after 
the center of gravity of the airplane has passed through the gust front. 
Therefore, since for most types of aerodynamic controls the center of 
resolution of the force will remain essentially fixed, the single compen- 
sating force in general will not he capable of canceling hoth the gust 
force and gust moment simultaneously. 

One such case may he said to exist, however; that is, when the gust 
force and moment huild up so rapidly as to he essentially steps. Exact 
cancellation is then theoretically possible if the control force is placed 
at the airplane’s aerodynamic center. This result (and also a result 
corresponding to eq. (23)) is in effect the solution given in reference 17. 
As indicated earlier in this report, the step approximation to the gust 
force and moment inputs should become increasingly valid as the flight 
Mach number increases sufficiently beyond unity, so that the analysis pre- 
sented in reference 17 may apply even more appropriately to flight at very 
high speeds than it does for the speed range considered therein. 


Application of Optimum Filter Theory 


Although the solutions given hy equations (23) and (25) are satis- 
factory from a theoretical standpoint, practically there are several 
objections to them. First, it is clear that the system’ s response must 
he hoth very rapid and precise in order to supply forces and moments 
matching those of the gust. Second, it may he argued that there is little 
point in attempting to huild a system that in effect cancels gusts of al 1 
frequencies, since, even if no alleviation is provided, the airplane’s 
inertia will prevent it from responding noticeably to gusts having fre- 
quencies very much larger than the airplane’s natural frequency. Since 
we presume to possess some knowledge of the probable distribution of fre- 
quencies in the gust velocity, an alternate approach is to try to lighten 
the task of the control system hy asking that it cancel only those fre- 
quency components of the gust velocity that actually do cause large normal 
accelerations and pitching velocities. The attractiveness of the Wiener 
theory of optimization is that it specifies just how to do this in order 
to satisfy a given criterion of excellence. 

Responses to arbitrary control- force and gust velocity inputs .- In 
the subsequent study, we again assume that a single control force is 
available for counteracting the forces and moments imposed hy gusts. 
Further, we shall assume that the step approximation to the gust func- 
tions is not valid so that the possibility of perfect alleviation is 
excluded. 
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Let us rewrite equations (2k) as 


sr(s) = Li(s)Ag(s) 
Q(s) = L 2 (s)Ag(s) 


v (26) 


where 


L x (s) = Gi(s)Z(s)+G 3 (s) 

L 2 ( s ) — H x (s) X(s)+H 3 ( s) 

It is clear from the form of equations (2 6 ) that the quantities L x (s) 
pinfl L 2 (e) may he interpreted as over- all system functions; they relate 
the outputs in 7 * (<?) and q(cp) to the gust input while the control is in 
operation. 

A linear relationship exists between L x (s) and L 2 (s) . By elimina- 
tion of X(s) it may he found to he 

L 2 (s) = p(s)L x (s) + ct(b) (27) 


where 


p(s) 

a(s) 


Hi(s) 

Gi(s) 

G l( s )B3 (s) -G3(s)H x (s) 
G x (s) 


The quantities p(s) and cr(s) are of course known, involving as they do 
the known system functions Gq(s), 3q(s). 

Now return to the physical domain. Let 


2 i(<p) = L [L x (s)] 
^ 2 ( 9 ) = L 1 [L 2 (s)] 
R(<P) = L _ 1 [p(s) ] 
S(<p) = L-^ffCs)] 


( 28 ) 



26 


NACA TN 3290 


From equation ( 27 ) 


^2(9) = L Ep( s )Li(s)+cr(s) ] 


bo that, letting 


P(<P) 


9 

= L" 1 [p(s)L 1 (s)] = r R(|)2 1 (9-l)4| 

Jq 


•we have 

Z 2 (9) = P(9) + s (9) 
Finally, from equations ( 26 ) 


7* (9) = J il(i)a-g(9-S)«i| 
-00 

l(9) = 2a(l)a g (q>-|)cL6 


(29) 


(30) 


(31) 


(32) 


It will he noted in equations ( 32 ) that we have used a more general form 
of the superposition integral than heretofore. The reason for this is 
that we wish to invoke the assumption usually made in random processes, 
namely, that the process has started infinitely distant in the past. 
Consequently, for an arbitrary origin of cp, the existence of inputs for 
negative arguments must he admitted, and this is accomplished in equa- 
tions ( 32 ) hy letting the upper limits of the integrals approach infinity. 

The error .- Equations ( 32 ) represent the variations in normal accel- 
eration and pitching velocity caused hy arbitrary inputs of gust velocity 
and the unknown compensating force. It is desired to minimize these 
variations according to some criterion when the gust inputs are random. 

We shall take as our error criterion a linear combination of the mean- 
square values of the variations 7 * (cp) and q(cp) over the infinite inter- 
val. We choose to minimize a combination of the motions in order to 
reduce the possibility of arriving only at solutions which minimize one 
motion at the risk of intolerably increasing the other. Thus, we seek 
to minimize the mean- square error 

^2 = T¥) + I (7 ,a ); a >0 (33) 

where a is a weighting factor, open to choice depending on which it is 
thought to be more Important to reduce, 7 ’ or q. 
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The choice of a mean-square error criterion is of course somewhat 
arbitrary. Criteria other than the mean square are possible and may lead 
to greater reductions in the variations 7 * (cp) and q(cp) . However, its 
use in connection with the gust alleviation problem is reasonable, in 
view of its property of weighting the minimization in favor of reducing 
the undesirable large errors while permitting many small ones. Also, 
its form is mathematically convenient, since it leads to an integral 
equation that can be solved by known techniques. 

Integral representation of the error .- Substitution of equations ( 32 ) 
in (33) gives for the mean- square error 


r r? -1 2 r 1 

= J 2a(i)ag( cp-£)d£ + i £ Zi(i)a g (q>-£)d£ 


Expanding the integrals, we have 


e 2 = 


1 

a 


m 


w 00 

J' 2 2 (lx) a g( c l>-£l)d|l J' Z 2 (| 2 )ctg(cp-£; 2 )d{; 2 J + 

OO 00 

[y ii(si)a g (cp-ii)diiy i 2i(^ 2 )n g (cp-i 2 )di 2 j (35) 


Now by definition 


f (cp) 


= lim 

$-^■00 




Then in equation (35) 


(36) 


- idha — f | Z 2 U 1 )d| i yz 2 U 2 )ag( 9 - ll )a g (^ 2 )d &2 + 

00 00 

~y Zx(li)<i|x J' Zi(| 2 )a. g (cp-| 1 )a g (cp-^ 2 )d£ 2 J dcp (37) 
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Tiie limiting process acts only on ag(<p-|j_) since all other terms are 
functions only of £i, | 2 . We define 

4> 

3dm i P a g(9"£i) a g(9 - £2)'3-9 = W(| 1 -| 2 ) (38) 

$-^•00 2 <t J 0 ° 

-<b 

This is the autocorrelation function of the random distribution of gust 
disturbance s. We assume that the function exists, and moreover, that it 
is known. Carrying the limit through the integrals in equation (37) and 
using (38) we get for. the mean-square error 

03 00 

6 s = F IzU x)*Si /'fe(i 2 )W(ix-la)d6 2 + 

'-'O V/Q 


00 00 

\J Zi(li)d|i^ > Zi(& 2 )W(i 1 -i 2 )di 2 


(39) 


Note that the mean-square error does not depend on the gust velocities 
themselves, hut rather only on their correlation function. 

Equation (39) can he cast in terms of one unkno wn function P(cp) 
hy means of equations ( 30 ) and (31)* Thus, in equation (30), let 


^ 0 ( 9 ) 



(md) 


Then 


9 

*x( 9) = P P(l)Bo(q>-|)a| (41) 

'-'o 

The upper limit cp is justified here since P(cp), R 0 (cp), and lx(cp) are 
impulse responses and hence must he zero for negative arguments. In 
addition, let 


U(!i) 



s(| 2 )w(ix-i 2 )ai 2 


Suhstituting equations (31) , (4l) , and (42) in (39) , we get 


( 42 ) 
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00 00 

= F s(6x)u(6x)d|i + 2 ( P(lx)u(6x)d| 

Jo Jo 


1 + 


00 00 

r p(ii)ciix r p(i 2 )w(ii-i 2 )d6 2 + 

wo «-/q 


£2 


/ <3-5 1 F P(o'i)Ro(^i-o'i)<3-cri j W(^x-^ 2 )d| 2 A* P(a 2 )R o (S2-cT2) < 3-0 r 2 

a Jo Jo Jo Jo 

m 


Now in the quadruple integral, reversing the order of integration results 
in 1 


00 00 CO 00 

J' P(ox) d cri J P ( 02 ) dcr 2 J' Ro(^i-cri)d|x J' Ro(l 2 -ff 2 )W(£x-£ 2) d £2 

a x cr 2 


Interchanging dummy variables (£x*= — 5*tTi, — s>tr 2 ) , we have 


00 00 00 00 

J P(gx)agx J P(52)<3-|2 J Ro(bx-Sx)<icrx J Ro(cf 2 “S 2 )W(ax-tT 2 )^cf£ 


£x 

00 00 


= F P(£x) d £x F P(^2)<3-52 F ^o(cTx)dcfx F Ro(ff2)W(£x-i2+ffx-ff2) d d 2 

Jo Jo Jo Jo 

(hh) 


It is clear from the second of the above forms. that the result of the 
double integration involving Ro(cp) is a function of (£x-£a) alone . 
Fi n ally, using the second of the forms (44), we combine the last two 
terms in equation (43) ; this gives 


00 00 
F P(£x)cL£x F 

vq 


“£ 2 ) + ~^Ro(^i) d ^i^ R o(^2)W(^i-| 2 +crx-o , 2)dcrJ 

(45) 


The bracketed quantity in (45) can be shown to possess all the properties 
of an autocorrelation function (ref. 24) . It is, in fact, necessary that 
this be true in order for our subsequent development to hold. Let this 
function be called ?(£x-£ 2 )* Then, returning to equation (43), we find 
that the mean-square error becomes 
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00 00 

* P s(| x )u( 6 i)a 6 i +2 f Pdijudjag 

v O vq 

co 00 

J P(li)dtx^ PdaMli-SaJaSa 


1 + 


m 


Integral equation.- Since S, U, and. Y are known functions (cf. 
eqs. (28), (42), and (45)) while P is unknown (eq. (30)), equation (46) 
is now of precisely the form considered "by Wiener in reference 18. As 
show therein "by an application of the calculus of variations, the mean- 
square error is a minimum if (and only if) P(cp) is constructed such that 
the following relation is satisfied: 

00 

u(q>) + f PdiMcp-i^dl! = 0 ; q> > 0 (47) 

'-'o 

Equation (47) is known as a Wiener -Hopf integral equation of the 
first kind for the unknown P(cp) . A general analytical method for solving 
it is described in reference 18 (cf., in particular. Appendix C of ref. 
18). Also, several numerical methods now exist (cf for example, refs. 

31 and 32 ) which can "be used to solve equation (47) approximately should 
an analytical approach prove impracticable. Hence, we can assume that it 
is possible to extract P(cp) from equation (47) • Having P(cp) , we may 
compute Z 2 (cp) from equation ( 31 ) . The Laplace transform of the control 
system necessary to achieve a minimum mean-square error is then obtainable 
from the relation (see eq. ( 26 )) 


_ L g (s)— Hs(s) 
%(s) 


m 


This is the sought-for quantity - the transfer function of the response 
in control force to gust inputs. 

The min im ized mean-square error, which results when the control- 
force system is constructed according to equation (48) may be found con- 
veniently by substituting equation (47) into (46) and using equation ( 31 ) • 
The result is 


00 


m in 



UdlHadlHSi 


(49) 
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Finally, in order to compare the minimized mean- square error with 
the error that results when no alleviation is provided (i.e. , X(s) =0 ) , 
it will he necessary to have the latter result. Setting X(s) equal to 
zero in equations (26) and using equation (39) > we have for the unmodified 
mean- square error 


where 


00 00 

= J hs(ii)d|i £ ha(^ £ Mti-^ S )dt 2 + 
00 00 

S3 ( £ 1 ) d-i 1 ^ S3(i2) w (ii - 52) (3 -la 


ha(cp) = L _1 [Ha (s) ] 
S3 (?) = L" 1 [Q3 ( s ) ] 


(50) 


Application of Results 


It remains to investigate some of the characteristics of a ‘control 
system that might he huilt according to the specifications of the theory 
just presented. To he acceptable, one should expect of such a system 
at least the following: (l) adequate stability j (2) practical feasibility; 

(3) a significant reduction in the mean-square error from the error that 
results when no alleviation is provided. Unfortunately, these require- 
ments are not necessarily compatible; in order to check them numerical 
analyses must be made, using specific airplanes and operating conditions. 

We have chosen as an example case a tailless triangular-wing airplane 
flying at a high subsonic Mach number. 

System functions.- Since we have determined in the first part of this 
paper that the stability derivative formulation is a valid appro xim ation 
to the exact equations of motion, we shall use it here to define the 
airplane system functions. Further, in the system functions G x (s) and 
Hi(s), we shall assume that the location of the control force remains 
essentially fixed. Thus, for a step control force IF, the control moment 
is -IFxq, where Xq is the (dimensionless) distance of the point of 
application of the force from the airplane center of gravity. With these 
assumptions, the system functions may be derived from equations (9) and 
(21). (The functions Gi(s) and Hi(s) are obtained simply by replacing 
c L (s) and 2^(3) in equations (9) respectively by l/s and -(xq/s) .) 

The results are 
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Gx(s) 


Gs( s ) 


s 2 r(s) 
" IF 


s 2 r(s) 

lOv, 


^s^e^+C^j+Xo ^OLg+pE^J " ^%,+Qj^ 


U* t(s+n) £ +k 2 ] 


[^ 2 - s (%+ c m^~ Qm a j+ C m g (s) |s (cl^Cl&J+ Cr^J j- 


%(s) = 


Bb(s) = 


where 


sQ(s) 

if' 


sQ(s) 

Ixi^ 


[(s+ti)^!^ 2 ] 

-s (Vxo+Cn^ - ^C^+XoC^ 


[(s4^) 2 4k 2 ] 


‘[cmgC^fn'a+C]^) - c Lg (s)(s<^ + <^J 
£t)' [(s+njZ+k 2 ] 


*' = ^ + 

1 ^ c I u ~ T l (^Qm^+Omd) + 

X/2 


U = 


<- 


25ti* 

n^+CL^-Cm^C^ 






>(5D 




It remains to define the gust functions Ci, g ( 9 ) and 0 ^( 9 ) whose 

Laplace transforms appear in G 3 (s) and B 3 (s) . In the absence of theo- 
retical results for these functions applicable to subsonic speeds, we 
shall assume they may be approximated by the expressions 


*\ 
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where p and r axe parameters. It is emphasized that these are not the 
gust functions of any known wing. They are considered to he physically 
reasonable, however, and, by virtue of their zero slope at cp = 0, are 
believed to be specifically representative of pointed-nose wings. They 
are strictly applicable only to flight at subsonic Mach numbers, in view 
of their asymptotic approach to steady-state values. The particular Mach 
number to which the functions correspond may be varied by changing Cr (co) 
Cn^C 00 ) } and r. As the flight Mach 
number is increased toward unity p 
should be reduced in value. The 
value of r should be chosen so that 
CjUgXcp) has its greatest positive 

value near the value of cp for which 
the airplane center of gravity passes 
through the gust front. The partic- 
ular functions which were used in 
numerical computations are plotted in 
sketch (d) . 


Finally, the Laplace transforms 
of equations (52), which are to be 
inserted in the expressions for 
G3 (s) and Ha(s), are 



Sketch (d) 


V s) - 

p 2 (s+p) -2rs 


> (53) 


°nig( s ) — ( on a ( 00 ) 


s(s+p) J 


Stability and the location of control force .- Having defined the 
system functions, we next consider the conditions which must be imposed 
on them to insure the stability of the airplane and control system. 

First, it is a necessity that the system functions themselves be stable. 
In order to insure their stability, it is sufficient to stipulate that 
their denominators contain no zeroes having positive real parts. This 
will be true in all four functions provided the damping factor p and 
the stiffness factor k 2 are positive. These requirements are of course 
the usual ones that arise in nearly all dynamic stability analyses. 


In addition, we must insure the stability of both p(s) and its 
reciprocal, since both enter into the analysis for the optimum control 
force (see eqs. (30) and (k0)). Hence, it is a requirement that both 
numerator and denom in ator of p(s) be free of zeroes having positive real 
parts. Now p(s) may be written in the form (cf. eqs. (27) and (51)) 
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_ H]_(s) _ 1 (s+cx) 

P S Gx(s) K (sS+^xS+^x 2 ) 

where 

K L_ 

Tl’xo+Cm^ 

^%t, +Cm a 

cx 

n’Xo+Cn^ 


^x “ £ 


( C m q +C i%) + *o ( G L a +G %) 


k x 2 = " 7 ( G m a +x o C I ta ) 


Then p(s) and l/p(s) are stable if c x , px> and k x 2 are all positive. 
We note in equation (54) that c x , Pi, and k x 2 are functions of Xq, the 
control location, and hence the requirements of stability restrict to a 
great extent the possible locations of the control force. Thus, from 
equation (54) 


kx 2 > 0 


if Xq < 


C La 


Pi > 0 


if Xq < 


c m q +C Pid. 

C ^ J q +C ^4d 


> (55) 


j 


With kx 2 > 0, we have 

Cx > 0 if Xq < - 

The condition k x 2 > 0 excludes as possible 
force al 1 points aft of the aerodynamic center x a ; whereas , with 
kx 2 > 0, in order to have c x > 0 we must exclude essentially al l points 
aft of the center of gravity. With k x 2 = c x = 0, the aerodynamic center 
is a possible location in the event that x a < - [ ( C^+C^) /( Cl^+C]^) ] . 

This cannot be considered a satisfactory position, however, since only a 
slight shif t of Xq in either direction from x a will cause cx to be 
negative in one case or kx 2 in the other. Hence, we conclude that in 


locations of the control 
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order to satisfy the requirements for stability, the control force must 
"be located forward of the center of gravity and sufficiently so to insure 
that xo < - ( Cja^J tj ’ ) • 

Gust correlation function .- One further piece of information needs 
to be supplied before the analysis for the optimum control system can be 
carried out; that is, the specification of the gust correlation function 
W(cp) . We shall use for this quantity an expression that lias been adopted 
by a number of authors in recent applications to gust analyses (refs. 10 
and 33 through 36 ) , namely 



Here, a.g 2 is the average intensity of the gust vertical velocities, Vg 2 , 
made dimensionless with respect to the airplane forward speed V. The 
quantity L is a measure of the so -ca l l ed ’’scale of turbulence” (ref. 34-) • 
Comparisons of equation (57) with the results of existing experimental 
studies have indicated that the value of L may vary from the order of 
several hundred feet to over a thousand feet. The reader is referred to 
references 10, 34-, and 36 for a more detailed explanation of the origin 
of equation (57) and the ranges of values to be expected of the quantities 

Wg and L. 

Numerical results .- Equations (51), (52 ) , and ( 57 ) complete the 
specification of quantities necessary for confutation of the optimum con- 
trol system from equation (-43) and the minimized mean- square error from 
equation (4-9) . Detailed results of the analysis are presented for the 
general case in Appendix B. In this section curves will be presented 
which illustrate the nature of the results as they apply to a single set 
of operating conditions. The numerical values which define the airplane’s 
inertial and geometric properties, its transfer functions, and the gust 
correlation function are listed at the end of Appendix B„ The set of 
constants defining the airplane are believed to be representative of those 
corresponding to a fighter-type triangular-wing airplane flying at a high 
subsonic speed. The constant L used in the gust correlation f uncti on 
corresponds to a scale of turbulence of 300 feet. The constant a,g 2 , 
defining the average intensity of turbulence, was left unspecified since 
it appears in the results for mean-square error only as a multiplying 
factor. Hence, the results can' be made applicable to a range of weather 
conditions by appropriate selections of this parameter. 

Mean-square error: Figure 7 is a comparison of the magnitudes of 

the minimized and unmodified mean-square errors (eqs. (4-9) and ( 50 )) for 
a range of values of weighting factor a and several values of force 
position. Large values of a correspond to weighting the minimization 
procedure in favor of reducing the mean-square pitching velocity, possibly 
at the risk of increasing the mean-square normal acceleration beyond its 
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■unmodified value. For small values of a the reverse is true. It is 
clear from inspection of the figure that the reductions in total error 
may he significant; the results are to he interpreted, however, as indi- 
cating only that a significant reduction in one or the other of the indi- 
vidual components of the mean-square error is possible. It remains to 
he seen whether such can he realized without the penalty of increasing 
one component heyond its unmodified value. This point is clarified in 
figure 8, which shows the individual components of the mean- square error 
given as ratios of their corresponding unalleviated values. (The note 
at the top of figure 8 permits one to attach absolute values to the indi- 
vidual components of the mean-square error, should this he desired.) It 
is clear that either error can he made arbitrarily small, hut eventually 
at the expense of increasing the other heyond its unmodified value. Of 
greater Importance is that a range of a does exist for which significant 
reductions in both errors are possible at the same time. Thus, for exam- 
ple, for the control force located at the wing nose (xq = - 0 . 25 ) , both 
errors may he reduced to about half their un all eviated values when a 
is near 0.1. 

In addition, note that the control force becomes even more effective 
in reducing both errors as its point of application is moved forward of 
the nose. This would suggest that a canard control surface might he an 
efficient means of generating the required force. It should he pointed 
out, however, that the analysis in its present form does not take account 
of the effects on the lift of the main wing of downwash from the canard, 
so that the analysis is applicable to the canard only if these effects 
can he considered negligible. As previously mentioned, the entire ques- 
tion of how to generate the required control force is considered to he 
heyond the scope of this paper, and the present results are of value 
mainly as an indication that more extensive analyses, in which this 
question must of course he considered, may he justified. 

Control force: We examine next the properties of the control force 

that achieves reductions of both errors, using as an example the case 
Xo = - 0 . 25 , that is, where the force is located at the wing nose. 

Figure 9 shows the transient response of the control force to a step 
gust. This force caused the mean-square values of normal acceleration 
and pitching velocity to he reduced to about half their unalleviated 
values. AIbo shown as dashed curves in figure 9 are the inputs in aero- 
dynamic force and moment due to a step gust. Note that the control force 
need he only about a third the magnitude of the gust force and need he 
applied much less rapidly. This should he contrasted with the problem 
of devising a system to provide perfect alleviation, for which the control 
force and moment must he made to match those of the gust at every instant. 

Further insight into the nature of the control force may he gained 
from inspection of the amplitudes of the control-force response to har- 
monic gusts. These are shown in figure 10. Also shown in figure 10 is 
a plot of the gust spectral density, normalized so that its value is 
unity at \ - 0 . (The spectral density is defined as the Fourier 
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transform of the gust correlation function, eq. (57) j it is a measure 
of the "energy” contained in each increment of frequency in the gust 
spectrum (cf. ref. 8).) Note that, as might he anticipated, the ampli- 
tude of control force falls off with increasing frequency as the gust 
spectral density does. Again, this should, he contrasted to the dashed 
curve which is the amplitude of control force required for perfect alle- 
viation. In the latter case the force is independent of the gust charac- 
teristics and clearly is of significant magnitude over a wider frequency 
range. 


Airplane responses: Figures 11 and 12 show the amplitudes of the 

responses in normal acceleration and pitching velocity to harmonic gusts 
with and without operation of the control. As in figures 9 and 10, the 
control force is located at the wing apex. On comparison of the results 
for a = 0.1 with the unmodified responses, it is seen in hoth figures 11 
and 12 that the amplitudes of the optimized responses have teen signifi- 
cantly reduced in the vicinity of the peaks of the curves and that the 
peaks themselves have teen shifted to the right where the gust spectral 
density (fig. 10) has begun to fall off. The effect of both modifications 
is of course to reduce the mean-square errors. Also shown in figure 11 
is the response in normal acceleration required to give zero pitching 
velocity response (a -> oa) and in figure 12, the response in pitching 
velocity which gives zero normal acceleration (a 0) . The reasons for 
the increase in mean-square error of one response, which is a consequence 
of attempting to minimize only the other response, are readily apparent 
from a comparison of the two curves with their respective unmodified 
responses. 

As a final point, it will be noted in figures 11 and 12 that the 
Initial values of the optimum harmonic responses are not zero. This 
means that the steady-state values of the corresponding transient responses 
to a step gust are not zero. Hence, in response to a step gust, the con- 
trol force will cause the airplane eventually to climb or dive. Unfortu- 
nately, the worsening of the transient response usually must be accepted 
as a consequence of attempting to mi ni mize the response to continuous 
random disturbances; one must count on the probability that isolated step 
gusts will not be encountered. 

Off -design performance: The results presented in figures 7 through 12 

apply only to a single atmospheric condition, namely, that corresponding to 
a scale of turbulence of 300 feet. It is of interest to ask how the system 
designed to operate optimally at one atmospheric condition will perform at 
others. To investigate this, some computations have been carried out for 
the control force which is optimum for a scale of t turbulence of 1000 feet. 
The results are shown in figure 13 in conparison with the optimum control 
force for L = 300 feet (from fig. 9) • It is apparent that, since the 
results are quite similar, the design is not strongly dependent on the 
magnitude of L. Hence, one may expect that a system designed to be opti- 
mum for one atmospheric condition will perform well for off -design 
conditions not too distant from the design condition. 
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The atmospheric conditions are also dependent on the quantity a.g 2 . 
However, this quantity occurs in the expressions for "both the minimized 

and unmodified errors as a scale factor. Hence, changing ctg 2 will have 
the same effect on "both of these errors and the errors of the op timi zed 
system will he reduced hy the same percentage for all values of a*g' 2 . 


Concluding Remarks 


The preceding application of the Wiener optimum filter theory to the 
problem of reducing the airplane’s response to random atmospheric turbu- 
lence has demonstrated that such an approach may lead to useful results. 

It was found that significant reductions of both the normal acceleration 
and pitching velocity responses may be achieved, and that construction 
of a control system, which realizes these improvements may be considerably 
simpler than construction of one designed to provide perfect alleviation. 
Nevertheless, the method and results presented here axe by no means meant 
to be final or definitive. For one thing, the problem for tailed aircraft 
has not been touched; the revisions and amendments of the analysis which 
are necessary to account for the tail and the effects of downwash may 
invalidate some of the results obtained here. Further, even for taillesB 
aircraft, in some cases the assumption that the control pitching moment 
is proportional to the control force may not be warranted. This would 
probably be the case, for example, in the event that the force is to be 
provided by a canard control surface whose interference effects are not 
negligible. Moreover, other formulations of the problem are possible. 

For example, rather than leave the control force unspecified as was done 
here, it may turn out to be practically more feasible to begin with a 
control system of known type having, say, open parameters and use the 
Wiener theory to give optimum values to these parameters. The main func- 
tion of the present work may be to offer evidence that efforts in these 
directions will likewise yield results that indicate the possibility of 
achieving significant reductions of the airplane’s response to random 
turbulence. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., June 28, 1957 
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■ APPENDIX A 

TRANSFORMATIONS TO HARMONIC. RESPONSES OF THE SUPERSONIC 
INDICIA! FUNCTIONS OF TWO-DIMENSIONAL, RECTANGULAR, 
AND WIDE TRIANGULAR WINGS 


The following is a compilation of the transformations to harmonic 
responses of the supersonic indicial lift and pitching-moment responses 
to step changes in (l) gust velocity, (2) angle of attack, and (3) pitch- 
ing velocity. The transforms are presented for the two-dimensional wing, 
the rectangular wing, and the wide triangular wing. Sources from which 
the indicial functions were obtained are listed as references 25 to 30. 


NOTATION AND TRANSFER OF AXES 


All coefficients to he presented are given as transforms of deficiency 
functions; these are defined as follows: 


Sx(iA) 

Ss(iA) 

fi(iA) 

fs(iA) 

f 2 (iA) 

f 4 (iA) 


00 

= J C Lg (o>)-C Lg (cp)J e" iA % 

00 

- Qmg(“)-Qmg(<p) J 

CO 

= f r^H-c^cp) 

w 0 L J 

co 

- J' c %' (°°) ~ c ma( ?) J 

00 

= J C L(i (oo)-CL (1 (cp)J 

00 

= J' C m ^(oo) -Qm^C cp) j 


-iAcp 
e dcp 


-iAcp 
e dcp 


-iAcp 
e dcp 


-iAcp 
e dcp 


-iAcp 
e dcp 




> (Al) 
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The transforms due to pitching velocity f 2 (iA) , f 4 (iA) will he presented 
for a wing pitching about the lea di ng edge or the apex. Likewise, all 
pitching-moment transforms are referred to an axis through the leading 
edge or the apex. To calculate the transforms for an axis location other 
than the nose, the following transfer relations may he used: 


Si(iA) = gio(iA) 

63 (iA) = gsoC^+iSxCiA) 

fi(iA) = fxo(iA) 

f s (iA) = f So (iA)+ifi(iA) 

f 2 (iA) f 2 o(iA)-|fi(iA) 

f 4 (±A) = f 4 o (iA)+if 2 o (iA)-5f 3 o (iA)-| 2 f 3.(11) 




} (A2) 




Here, | is the number of chord lengths of the new axis from the nose, 
measured positive rearward from the nose; the subscripted terms are those 
listed below, referred to the axis | = 0. 

TWO-DIMENSIONAL WING 
Transforms 


*(«) - 1 

[f 0 (M,to) -1] 

(A3) 

Bs ° (iX) 

[l-2f x(M,u) ] 

(A4) 

= A 

|x[f x (M,u) -f 0 (M,u) ]+i[f 0 (M,u) -l]j- 

(A5) 

f3o (iA) . | 

-j^A[fo(M,cj) -f 2 (M,w) ]+i[l-2f x(M,u) ]| 

(A6) 

i 

jx[2fi(M,m) -f 2 (M,w) -f 0 (M,u) ]-i[l+2f i(M,co) -2f 0 (M,u) ] 



(AT) 
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f4 ° (iX) = m{ Xf2(M ^ + GS " A ) ]+i[i-f 2 (M,^)+ 

3fx(M,u)-2f 0 (M,13)]| (Ag) 

In the above, the f n (M,t3) are Schwarz functions (ref. 37), defined as 


x 

k (M,w) = r u n e" iuu J 0 (^) du 
'-’o V 11 / 


- _ 2 &l 
w ~ P 2 


(A9) 


Numerical tabulations of these functions for M ^ 1.20 may be found in 
reference 38 . ' 


Initial Values 


As an aid in plotting the transforms as A 0, we give below their 
first-order expansions. 


_ gi(iA) 

A“»0 


2M 2 

P 3 


irn 

3P 3 


(2M 2 +l) 


(A10) 


SsoCiA) 


- 5 $ + T§ (a**) 

3 pa 4pS 


i^o fx(iV) 


_2_ . fw 
P 3 " P 3 


(All) 

(A12) 


lim f3o(iA) 

A-^0 


_4_ + 3 iu 
3P 3 4 p 3 


(A13) 


^S 0 W») - 


2 

3P 3 


iui 

4P 3 


(A14) 
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m 0 ^o(») 


1 ±03 

23® + 53 s 


(A15) 


Stability Derivatives 


Finally, in order to complete the compilation of terms needed to 
calculate the harmonic responses, equations (13) and (l4) in the text, 
the necessary stab ilit y derivatives are listed helcrw. 


%,(") = ( V 00) = 1 

(Al6) 

Qn k 0 (oo) = Pm^C") = - | 

(AIT) 

§ 

11 

1 

(A18) 

/ 

11 

U) I 

(A19) 

< W“ ) * I 

(A20) 

s>> - - 55 

(A21) 


RECTANGULAR l-JXNG 


The harmonic functions and stability derivatives for the rectangular 
wing are each composed of a two-dimensional, contribution plus a contribu- 
tion dependent on aspect ratio. Listed below are the aspect-ratio contri- 
butions, each of which should be added to the corresponding two-d im ensional 
function, already given. The results are valid only for those combinations 
of Mach number and aspect ratio for which f3A ^ 1. 


Transforms 


-SL ( 2 + + 

MAA S \ M ^2 J 


dgx(iA) = 


(A22) 
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Ag3 D (iA) 

Afi(iA) 

Afso(iA) 

Af 2o (iA) 

Af 4o (iA) 

where 


MAA 3 


A + it i E? + 2 ±( me-iV^-iAM 
LMA 3 V J 


| ^ e -iA/m_ n2e -iA/n^ 


(A23) 


m 2 aa 2 


fc^AA 2 


-2 + i ^2^ - 0 i ^me“^/ m +ne“^^/ n ^ (A24) 

1 + ~ i ^ i ^m 2 e“ i A/ m +n 2 e - i A/n\ _ 

L M 2 ^ 2 33 2 A 2 \ J 


± ^e-^+ne’ 1 ^ 


(A25) 


l^AA 2 


- 1 + + i (p . |) . £ (m-e-»^e-W=)' 


(A26) 


3m 2 aa 2 


2 + 4i 


1 - If2f + 2(M 2 +3) 


4p £ 


M^ 2 


+ ~ ^m 2 e “ 1 ^/ m +n 2 e " 1 ^/ n ^ 


^#e _i V m +n 3 e -iA/n^|. 


(A27) 


n 


Mt-1 

M 


m 


M-l 

M 


Initial Values 

+ gpi (3M 2 + 1) (A28) 
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11m A&s 0 (iA) 

A— ^0 


if. - 1. US. (3M 2 +l) 

3 4 a 15 p 4 A 


Urn Afi(iA) 
A-> 0 

11m Afs 0 (l-A) 
A->0 


11m Ztf2o(lA) 
A-»0 ^ 


11m Af4 n (lA) 
A -^-0 0 


2 (1^+1) + ~^ 7 ~ (M 2 +3) 


3P 4 a 


6p 4 A 


— i- (M 2 +l) - — (M 2 +3) 

2f3 4 A 15 P 4 A 

— L_ (1^+1) + -i- (m 2 -^ 

6p 4 A 30 p 4 A 


_2 (M 2 *!) _ 2 l (^+3) 
15 f3 4 A 36 P 4 A 


(A29) 

(A30) 

(A31) 

(A32) 

(A33) 


Stability Derivatives • 


ACiJco) = ^(00) = - ^ 

(A34) 

- ^go (oo) - 3 p 2 A 

(A35) 

ACitf, = (^+1) 

Lx 3P 4 a 

(A36) 

ACnido = " ( M2+1 ) 

23 4 A 

(A3T) 

A( W” ) ’ ' 3? 2 A 

(A38) 

% w = 2 P 3A 

(A39) 


■WIDE TRIANGULAR WING 


Presented below are the transforms, their initial values, and perti- 
nent stability derivatives for the wide triangular wing. Results are 
applicable so long as the wing leading edge is supersonic; hence it is a 
requirement that 3 A > 
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Sx(iA) 

&3 0 (iA) 

f“x(iA) 

^oCiA) 

fsodA) 

f4 0 (iX) 


Transforms 


±[2 

PA (A 


(i, g 


-i 


8 

3PA 


{*[ 


5-i ^o ( 5 , jp) -fl(M,S) 




m 


= ^ {a 2fi(M,(3) -f a (M ,to) -f 0 (M,w) 


+ 


2f 0 (M,T3) -2f oXM/ui ) -1 


} 


= 3PA + ( X “•£) f o(M,u5) -2f x (M^ 


l-2f o (M,0))+3fx(M> w ) “^ 2 (^ 10 ) 




-l-2f 2 (M,u)+2f 0 (M,w) 


2f 1 (M,u)-f 0 (M,u) 


(a4q) 


(Aia) 


(A42) 


(A43) 


(A44) 


: A + ^i) -2 (a + I) fx(M,75) + 

f pm 2 V\ r 1 

i ^ 1 + y Af 0 (M,ui) + i|l+2f a (M,fi) + A fi(M,u) 

2 ( 1 + s) 


(A45) 
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The f Q Qjj, function which appears in 
and (A4l)) is defined as. 


gi(iA) and gs(iA) (eq.s. (A*K)) 


fo 




(A46) 


Tabulations of this function are not as yet available; in their absence 
the function can be evaluated either by numerical integration or by means 
of the series expansion, (cf. ref. 37) 


(* w) = e 112 1 (• + 1Jn+1 

n=o ■ , . 


(A47) 


Initial Values 


#So *<») ■ £ - g? «« +s) 


u-, •<.(«) - - £ 


(A lt 8 ) 
(A49) 


iSo fl(ix) ' sF ■ # 


/x 1 2 ito 

7 JS 0 f3o(lA) = V + 5 p 3 


aSo f2 ° (iA) " 3^ " 5P 3 


(A50) 

(A51) 

(A52) 
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Stability Derivatives 


%.<”> = °LgW - J 

(A54) 

°^ ( ") = < W“> * - •§ 

(A55) 

n 4 

14 ' 3B 3 

(A5 6 ) 

/ 

II 

^l‘H 

(A57) 

Cl Qo^“ I 3.8 

(A58) 

So'" 1 = ' I 

(A59) 
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APPENDIX B 


SOLUTION OF INTEGRAL EQUATION FOR P(cp) 


The following is a compilation of expressions -which arise in solving 
the integral equation for’ P(cp) , equation (47) of the text. The equation 
is readily solved hy the use of the two-sided Laplace transform (ref. 39) • 
For any function F(cp), define the transfqrm as 


00 

= J F ^ C P) e ” BCPdc P 
-00 


(Bl) 


and we denote the transforms that result from integration of F(cp) for 
positive and negative values of cp hy f + (s) and f_(s), respectively. 
Thus 


f + (s) = 

00 

r e“ s ?F(cp)dcp 

V 

Jo 


0 

f-(s) = 

J ' e*’ s %(cp)dcp 


-OO 




> (B2) 

- 


The solution may he put in the form 


- / a({p)e' BC Pdcp 

p+(s) = ~ 2 l&) (B3) 

where p + (s) is the transform of the desired, quantity P(cp) • The latter 
quantity is different from zero only in the range cp > 0; it is zero for 
cp < 0, as it must he to he physically realizable as an impulse response. 

The following quantities arise in the derivation of a(cp) and 0 + (s) : 

cr + (s), S (cp) : From equations ( 27 ), (51 ) , and (53) 


Ags^-AxB 


(B4) 


a+(s) = 


(s+a 1 )(s+a 1 *)(s+p) 3 



where 


The inversion of 
S(cp) 

where 


ax “ M-i + i ’^ 1 2 -|A 1 2 





(^x +X ° G 0 


•2rCma 


Ax = ^- (^C^+XqCi^ 
cr + (s) gives 


= (so+SiCfH-Saq) 2 ) e"*^ + S3e“ ai ^ + s 3 *e -a:L * CI) 


s 2 = 


sx = 


7vgS 2 +Ax^ 

2 L(s+a x ) ( s+ai*) J g = _ p 


ds 


A2S 2 +Ax s 


L(s+a 1 )(s+a 1 *) J g = -p 


1 d 2 
s o = r ~r ~2 


AaS^Axs 


s 3 = 


S3* = 


2 ds 2 L(s+a 1 )(s+a 1 *)J s=: _ p 
AgS^+AxS 


L(s+p) 3 (s+a 1 *) J s _ _ 




ax 


Aa s2 +^x s 

L(s+p) s (s+a 1 ) J s = _ ai * 
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p + (s) : From equations ( 27 ) and ( 51 ) , 


P+(s) 


1 s+Cx 

K ( B+&x) ( s+ai*) 


•where 

K - - i 

n'Xo+Cny. 

x o G rtt +Cm a 

Cl = 

n’^o+Cn^ 

w(s) : Using equation (57) 


■where 


w(s) 



U(cp) e“ sc Pdcp 


= 3Na g 2 


y 2 - s 2 

'(I'^-s 2 ) 2 


(b6) 


(B7) 


N 


e 


L 


7 s 


3 


u(b) , U(cp) : Using equation (42) 

00 

u(s) = J' U(cp) e“ sc ^dcp = cr + (s)w(s) 


(B8) 


To find U(cp) for cp > 0, we evaluate the residues of the product of 
cr + (s) and w(s) at poles in the left -half plane. The result is 


-pep 


+ U 3 e 


-a x cp 


+ 


U(cp) (cp > 0) = (Uo+UicpfUacp 2 ) e 
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where 


U 0 = s 2 w**(-p) + 8iW*(-p) + s Q w(-p) 

Ux = 2s 2 w*(-p) + Siw(-p) 

U 2 = s 2 w(-p) 

U 3 = S3W(-a!) 

U 3 * = ssM-ao.*) 


U 4 * a g 2 


ff+(-K) - | tr + ' (-N) 


U 5 =-<x g 2 | ff + (-H) 

i|/-(s): The two-sided Laplace transform of the autocorrelation function 

Y(cp) may he -written (using eq_. (h5)) 

00 

i^-(s) = P ¥(cp)e -s< ^dq) = w(s) 1 + — — 

L ap.(s)( 

K^s) 


a P + (s)p + (-s)J 


a( Cl 2 -s 2 ) 


( s 4 -27 2 s 2 +y 0 2 ) 


(BIO) 


where 


2y z = a x 2 + ax* 2 + •— 
7 0 * - ( a ^*) 2 + 

This may he factored to the desired form 


1 


i}r(s) = e + (s)e_(s) 


(Bll) 
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where 


( 7 +s) (h+s) (h*+s) 
(K+s) 2 (c 1 +b) 


9 ( s ) - g <£! 1 (b - e) (b *- B) 
" (N-s) 2 (ci-s) 

and 



a s 


a(cp) : Letting 


A(s) 


^(s) 

9_(b) 



(B12) 


we find a ( 9 ) for cp > 0 “by evaluating the residues of A(s) at poles in 
the left-half plane. Let 


A_(s) 


1 

0-(s) 


(B13) 


Then 


a(cp) (9 > 0) = (A 0 +A 1 qH-A 2 9 2 )e PC? + A a e~ a ' 1 ^ + 
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•where 


P(<p) ; 

find 


where 


Ao 

= UoA^C-p) + 

UxA_ * 

(-P) 

Ax 

= UiA^-p) + 

2UaA. 

.*(-p) 

As 

= U 2 A_(-p) 



A3 

= U 3 A_(-ax) 



A 3 * 

= U 3 *A_ ( -a^*) 



A4 

= u 4 Aj-N) + 

u 5 a_* 

(-N) 

A5 

= u 5 a>n) 




Having determined 0 + (s) (eq.. (Bll)) and cx(cp) (eq.. (Blh)), we 
P(cp) hy inversion of equation (B3) . The results may he written 

P(<P) = ( p o+ I> i ( l } i I, 2 ( P 2 ) e *" PCp + p 3 e "" ai< ^ + P 3 *e“ ai * cp + 

Pae-” 1 + P 7 e- bt P + 


Po = -UoA(-p) - UxA'C-p) - U 2 A l, (-p) 

P x = -UjAC-p) - 2UaA*(-p) 

P 2 = -UaA(-p) 

Ps = -H3A(-ax) 

P 3 * = -U a *A(-a 1 *) 


p e = 


(H-y) 2 (ex-y) 
(h-y) (h*~y) 


A + (-7) 


P 7 = 


1 ( _ t) 

* (7-tOCb*-*) 


P 


* 

7 


1 (N-h*) 2 (c x -h*) 
H (y-h*) (h-h*) 


Af(-h*) 
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A(s) 


A + (s) 


1 


t(s) 



io + 

2 Ae 

^ + 

Pts 

(p+b ) 2 (p+b ) 3 

(ai+s) 


As* 

+ J± + 


a x *+s 

N+s 

(N+s) 2 


TABULATION OF CONSTANTS 


Given Below is a tabulation of constants which were used in numeri- 
cal computations to describe the airplane’s inertial anfl geometric 
properties, the system functions, and the gust correlation function. 

Airplane : 


Mass parameter t)’ 
Inertia parameter £ 
Wing chord c 
Center of gravity located 


= 100 
= 100 
= 30 ft 
at 0.25c 


Airplane system functions; 



3.0 

Si 

= 0.50 

G ns a = 

-0.75 

H 

= 0.0175 


2.0 

k 

= 0.0861 

= 

-1.0 

P 

= 0.20 


-1.0 

r 

= 0.021 


Gust . correlation function: 


L = 300 ft 
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Figure 1.- Normal acceleration response to harmonic gusts of an aspect- 
ratio-3 rectangular wing flying at Mach number 1.2. 
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Figure 3 .- Normal acceleration response to sharp-edge gust of an aspect 
ratio-3 rectangular wing flying at Mach number 1.2. 
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Figure 4.- Pitching velocity response to sharp-edge gust of an aspect- 
ratio-3 rectangular wing flying at Mach number 1.2. 
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(a) B.P. [Zl (») 



(b) I.P. ZL (i A) 

L a s 


Figure 5 .- Normal acceleration response to harmonic gusts of an aspect- 
ratio-3 rectangular wing that is restrained from pitch i ng; Mach 
number 1.2. 
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Figure 8.- Ratios of Individual components of minimize 

components of unmodified mee 
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Figure 10.- Amplitude of control-force response to harmonic gusts. 
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Figure 11.- Amplitude of normal acceleration response to harmonic gusts 
with and without operation of control force. 














